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Abstract. In April 2001 SN1993J was observed with both the PN and MOS cameras of the XMM-Newton 
observatory, resulting in about 7 x 10^ s of acceptable observation time. Fit results with both the PN and MOS2 
camera spectra studying different spectral models are presented. The spectra are best fitted in the energy range 
between 0.3 and ~ 10 keV by a 2-component thermal model with temperatures of kTi — 0.34±0.04 keV and 
kT2 — 6.54±4 keV, adopting ionization equilibrium. A fit with a shock model also provides acceptable results. 
Combining the XMM-Newton data with former X-ray observations of the supernova, we discuss the general trend 
of Lx oc t~'''30 and the bump of the X-ray light curve as well as former and recent spectral results in the light of 
the standard SN model as first proposed by Chevalier in 1982. 

Key words. X-rays - supernova - SN1993J 



1. Introduction 

M81 at the 



(Freedman et al. 1994) 



moderate 
is 
of 



distance 
one of the 
March 1993 



of 3.6 Mpc 
best observed 
the supernova 



galaxies. When end 
SN1993J exploded in that galaxy, there was good hope to 
learn about the progenitor and follow, as with SN1987A in 
the LMC, the evolution of the early phases for years with 
high quality spectra in many wave bands. 5 days after the 
detection in the optical bands ( (Ripero & Garcia 1993) I, 
radio emission from SN 
GHz 



(Weiler et al. 1993) 



1993J was observed at 22.4 
by the VLA, demonstrating 
the interaction with circumstellar material (CMS), an 
important prerequisite for X-ray emission. 1 day later, at 
day 6, the SN was observed with the ROSAT satellite at 
soft X-ray energies around 1 keV (12 Angstrom). The im- 
mediate detection of a strong X-ray signal (Zimmermann 
et al. 1993a, b) at these early times was unprecedented 
up to that date, even though in the following years 
a few of the 18 presently known X-ray supernovae^ 
have been detected at even earlier epochs. The X-ray 
satellite ASCA with an extended spectral coverage 
(Tanaka et al. 1993) and the hard X-ray instrument 
OSSE onboard GRO [(Leising et al. 1994)] extended a 



few days later the spectral bands in which emission of 
SN1993J was observed. 

X-ray observations of the supernova were continued 
with the ROSAT observatory up to 1998 and with the 
ASCA satellite up to 1996. After a break of 2 years the 



CHANDRA observatory observed SN1993J in May 2000 



(Swartz et al. 2003)| followed by the XMM-Newton obser- 
vation in April 2001, reported here. 



2. The Supernova Scenario 

The fact that 93J began as a type II SN, characterized by 
prominent H lines in the optical spectra and few months 
later changed to the lb type, characterized by no or few 
hydrogen but rich in He, is interpreted that the progenitor 
had lost all but a small amount of its hydrogen layers, thus 
that the shock heated effective photosphere could quickly 
sink through the thin H layer into the deeper He layers 
during the initial expansion and cooling phase. 

In the case of SN1993J, optical plate material taken 
before the explosion, the optical light curve and the- 
oretical considerations indicated, that the progenitor 
was probably a KG I red supergiant with an ini- 
tial main sequence mass of about 15 Mq , that had 
kept less than 1 Mq of H at the time of the explo- 
sion 



(Nomoto et al. 1993) (Podsiadlowsky et al. 1993) 



(Ray et al. 1993)||(Bartunov et al. 1994),.(Utrobin 1994) 



Send offprint requests to: H.-U. Zimmermann 
^ http:/ /www. astro. psu.edu/~immler/supernovae_hst. html 



(Woosley et al. 1994)| [(Aldering et al. 1994) 

Loss of the outer hydrogen layers of the star may have 
been caused by either a strong stellar wind or by mass 
(Roche lobe) overflow to a companion star in a binary 
system. It is known that red supergiants blow slow winds, 
in the order of 10 km/s and carrying between 10^^ to 
10^** Mq per year, into space. It is unclear whether a 
wind alone can account for the extreme loss of hydro- 
gen in a star of 15 Mq, therefore there is some tendency 
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to prefer a binary scenario for SN1993J with mass over- 
flow to a companion, while additional material escaped as 
a wind. In a binary scenario one might expect disk-like 
asymmetries, which indeed have been found in the pro- 



files of some optical lines (Matheson et al. 2000) On the 
other hand radio images of the SN1993J explosion neb- 
ula (Bietenholz et al. 2001) show nearly perfect circular- 
ity with no indications for asymetries (< 6%). 

In the coUaps of the inner core of a supergiant star 
an extreme Shockwave traverses the layers of the super- 
giant, driving the upper layers of the star into a circum- 
stellar material producing strong shocks in that medium. 
Observations of the shock interaction in the radio, optical 
and X-ray bands provide inputs for understanding the ex- 
plosion dynamics as well as element composition and mix- 
ing. Its physical basis is usually described in the standard 



or 2-shock model first worked out by (Chevalier 1982) 



In this model the expanding supernova drives a for- 
ward shock with a typical velocity of 10"* km/s into the 
circumstellar matter giving rise to temperatures in the 
order of 10^ K, ~ 100 keV. The circumstellar matter is 
assumed to originate from the stellar wind of the super- 
giant progenitor star. The resultant density profile is then 
a function of the mass loss rate by the wind, the wind 
velocity and has a radial dependance of r~*, where s=2 
if the wind is isotropic and constant. In the model also a 
reverse shock is forming by the interaction of the shocked 
material with the SN ejecta, from which temperatures in 
the order of 1 keV are expected. Initially the reverse shock 
region may be radiative and can form cool material that 
may absorb the X-radiation of that region. The density 
profile of the SN ejecta is also described by an r~" depen- 
dance. One can then derive a self similarity solution for 
that scenario where n and s, describing the density pro- 
files of the ejecta and the circumstellar medium, play an 
important role. 

The reason that early X-rays are so rarely seen, up to 
now only from 18 events, is mainly caused by the circum- 
stellar environment. If the density is very high, then the 
initial radiation may not be able to ionize all the circum- 
stellar material and X-rays may be fully absorbed. If, on 
the other hand, the circumstellar matter density near to 
the SN is very small then the intensity of the X-rays from 
the shocks may be too weak to be detected with present 
day instruments, as was the case with SN 1987A. 

3. XMM-Newton observations 

XMM-Newton was pointed towards the center of M81 - 
SN1993J is at 3 arcmin distance - on April 22/23, 2001, 
for a total of 132 ksec, from that data of 90 ks with the PN 
camera in small window mode and 83 ks with the M0S2 
camera in image mode were taken. The medium filter was 
applied both in the PN and the M0S2 instrument setup 
during the observation. The MOSl camera was set to tim- 
ing mode and not used in this analysis. 

Due to high and variable particle background in parts 
of the observing period, a temporal screening had to be ap- 



plied. While for the M0S2 camera only the last few 1000 s 
were screened off, leaving 77 ksec of good data, in the PN 
observation about half of the observation was affected by 
bursts in the particle background. In order to maximize 
the signal to noise ratio in the PN spectrum, we deter- 
mined the rates in intervals of about 1000 s separately 
for the source plus background field and the background 
field and eventually excluded all time intervals that de- 
teriorated the signal to noise ratio in that band, leaving 
69 ks of good data. The optimization relative to the cut 
radius of the source plus background field and the deselec- 
tion of background affected time intervals was performed 
in different energy bands. The final selection of accepted 
intervals was done in the band above 1 keV to optimize 
the signal in the part of the spectrum where the lines 
clearly stand out in the spectrum. The cut radius chosen 
was 300 pixels, corresponding to 15 arcsec. The spectrum 
was binned such that each bin contained at least a fixed 
number of 5 (background subtracted) source counts. 

4. The XMM-Newton spectrum 

The PN spectrum of SN1993J shows emission lines of 
highly ionized Mg, Si, S, Ar, Ca and very cleary the com- 
plex of the Fe hnes at high energies. Spectral fits were first 
performed using only the PN camera data offering the best 
statistics. Fitting was tested with different model com- 
ponents. Good fits were achieved over the whole energy 
band between 0.3 up to ^ 11 keV, using a 2-component 
thermal model with variable element abundances (vmekal 
in the nomenclature of XSPEC). We also tried different 
shock models and found that a sedov model with variable 
abundances (vsedov in XSPEC) produces also acceptable 

values, although all tested shock models tend to repro- 
duce not so well the different line complexes visible in the 
spectra above 1 keV. 

An overview on the best fit parameters with both the 
PN and the M0S2 camera and for different models is given 
in Tab. ^ The PN camera spectrum and the best fit 2- 
component thermal model (2vmekal) is shown in Fig. ^ 

The spectral differences found between the PN and 
M0S2 fits can be understood as uncertainties in the spec- 
tral cross calibration of the two instruments (as stated 
in the document 'XMM-EPIC status of calibration and 
data analysis', XMM-SOC-CAL-TN-0018 2003, valid for 
the software version SAS 5.4 applied in this analysis). In 
the following we therefore take the best fit results from 
the PN camera with its better statistics (7515 counts com- 
pared to 2593 from M0S2), especially at higher energies 
more representative for the overall behaviour, as the basis 
for the spectral discussion. 

The Fe complex near to 6.7 keV is represented in the 
PN camera with about 33 counts (in the MOS camera 
about 5 counts). The width of the Fe complex in the PN 
spectrum (extending from 6.5 to 7.2 keV) is generally bet- 
ter fit with a higher 'high temperature' component, while 
other line complexes, or better the regions between them 
tend to be better fit by somewhat lower temperatures. But 
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Table 1. Selected best fit parameters of different models. 2vmekal is a 2-component thermal model with variable 
element abundances, 2vmekall is the same model but with identical element abundances for the 2 temperature com- 
ponents, and vsedov is a shock model with variable abundances. 



instrument 


model 


kTi[keV] 


kT2[keV] 


iVffi[10^^cm"^] 


NH2[W'^^cfn'^] 


dof 


2 

Xr 


PN 


2vmekal 


U.Oy:_o.03 


6.54^1 


0.55±0.21 


0.33±1.80 


360 


0.90 


PN 


2vmekall 


"•'JJ_0.04 


4.241?:* 


0.45±0.13 


0.18±0.75 


371 


0.94 


PN 


vsedov 


0.51 ±0.04 


7.0 ±0.2 


0.17±0.037 




376 


1.01 


M0S2 


2vmekal 


50+0 03 

"•"-'^-0.04 


3.48^1:^ 


0.25±0.25 


0.0±0.54 


113 


1.25 


PN+MOS2 


2vmekal 


r, qq + 0.46 

"•■J'J-0.46 




0.59±0.25 


0.2±0.8 


505 


0.99 



10 



channel energy (keV) 



Fig. 1. XMM-Newton spectrum as observed with the PN 
camera and the best fit thermal 2-component model with 
fcTi=0.34 keV and kT2=6.5A keV. The dashed curves show 
the low and high temperature component. 



Table 2. Best fit element abundance values from the fit 
to the PN camera data 



cl. 


2vmckal_Tl 


2vmckal_T2 


2vmokall 


MOS2.T1 


vsodov 


Ho 


15±95 


0.1±273 


0.3±12 


0.0±156 


259±4881 


C 


0.0±172 




O.OiO.l 


0.0±156 


0.0± 


N 


0.0±172 




140±553 


212±1908 


6.5±120 





0.15±1.1 




0.2±0.9 


O.OilO 


1.1±20 


No 


2.0±11 


0.0± 


2.5±10 


0.0±14 


4.0±75 


Mg 


0.4±2.6 


47±4531 


0.2±1.0 


2.0±16 


6.4±119 


Al 


0.0± 


103±10000 


0.0± 


13±104 




Si 


3.6±20 


20±1930 


4.7±20 


3.7±30 


7.0±131 


S 


15±85 


10±968 


8.4±35 


8.2±68 


0.5±10 


Ar 


132±736 


18±1777 


16±68 


0.6±43 


5.8± 


Ca 


95±541 


3.7±357 


11±46 


18±169 


0.0±21 


Fc 


1.1±5.9 


1.3±121 


1.0±4.0 


1.8±15 


4.8±90 



fit requires additional intrinsic absorption. The lines at 
higher energies are not very well reproduced by the model. 

5. Discussion 

How do the X-ray observations fit into the scenario for 
SN1993J ? 



the statistics in the spectrum do not reasonably suggest to 
introduce more components with additional temperatures. 

Let us first look at the results from the 2-component 
thermal model fit (2vmekal). The low temperature com- 
ponent of kTi ~ 0.34 ± 0.04 kcV dominates the spectrum 
up to energies of about 2 keV, but a higher temperature 
component is definitely needed to explain the spectrum at 
higher energies. 

Given the fact that the high temperature component 
takes into account just the band between 3 keV and ^^8 
keV, where the spectral fit anchors on the Ca IXX, XX and 
the prominent Fe XXIII, XXIV, XXV line complex, the 
temperature is likely to be not very well defined. Indeed 
the temperature of the high temperature component of 
kT2 ~ 6.4 ± 4 keV is poorly restricted, because every 
change in the temperature component is compensated by 
changes in the abundances of Fe, Ca and Ar. Due to 
this dependence all the element abundances in the fits 
have errors typically much larger than the parameter val- 
ues themselves. Therefore, the abundance values should be 
considered with some care. 

One major difficulty with fitting the vsedov model is 
based in its enormous demands on computational power. 
It is therefore difficult to properly scan the parameter 
space. In the vsedov model the best fit temperatures are 
somewhat higher than those in the vmekal models. The 



5.1. Earlier X-ray spectra 

We have reevaluated the ROSAT data in a homogeneous 
way'(Zimmermann et al. 1997) and fitted the spectra with 



a thermal 1-component model (vmekal in XSPEC) where 
the element abundances were fixed to either solar or to 
the values obtained from a similar fit to the XMM-Newton 
PN spectrum. Due to the high temperatures in the early 
ROSAT observations there is almost no difference between 
fitting either solar or XMM-Newton based element abun- 
dances. In the November 1993 measurements the fit with 
the thermal model and solar abundances turned out to 
be poor but improves significantly (see Tab. |3Jl with the 
XMM-Newton based element abundances. 

The first ROSAT PSPC spectra do only allow to set 
a lower limit of > 17 keV for the temperature of a ther- 
mal brcmsstrahlung, Raymond & Smith or mekal model. 
From the ASCA measurements a lower limit of > 10 keV 
was reported. Around the same time the hard X-ray in- 
strument OSSE on the GRO satellite could determine a 
temperature of 82^29 keV at day 10 to 15 and, statisti- 
cally already very weak, of 7414q" keV about 1 month af- 



ter the event (Leising et al. 1994), ROSAT measurements 
half a year later (Zimmermann et al. 1993d, 1994a) re- 
vealed a strong decrease from the initial high tempera- 
tures to temperatures around 1 keV. An ASCA measure- 
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Table 3. ROSAT spectra: best fit parameters of the fit 
to a l-component thermal model (vmekal in XSPEC). 
Abundances were fixed to either solar or to those that re- 
sulted from a fit to the XMM-Newton PN spectrum using 
the same model. 



date 



April 93 
April 93 
Nov. 93 
Nov. 93 
April 94 
April 94 



kT [keV] NhIW^ 



> 17 

> 18 
1.601°-, 

3+0.5 
-0.3 

1 05+°-^^ 



0.56^ 



U.U/0_Q 

071+0 005 

J^-0.015 
"•UO' -0.015 

14+003 

"•^^-0.06 
'^-0.03 

n 14+003 



'] dof x2r 



abund. 



0.97 
0.96 
2.46 
0.90 
1.05 
0.90 



solar 

PN 

solar 

PN 

solar 

PN 



ment (Uno et al. 2002) confirmed that strong temperature 
drop. 

In terms of the standard 2-shock model it is suggested 
that the observed emission initially originated from the 
fast forward shock, while emission from the reverse shock 
region was blocked at that time by absorbing material due 
to fast cooling processes in the denser environment. It has 
been assumed that the initial absorption disappeared on a 
time scale of the order of 100 days so that half a year later 
the measured flux was then dominated by the emission 



from the reverse shock region ^(Zimmermann et al. 1994) 



(Zimmermann et al. 1996)1 (Fransson et al. 1996) 

Unfortunately the X-ray lightcurve of SN1993J is not 
covered between days 50 and 200, so that the transition 
from the forward shock dominance to the reverse shock 
remains an open issue. 

The Chandra spectrum of May 2000 was fitted by 
(Swartz et al. 2003) with a 3-component 2vmekal-|-mekal 
model. Their best fit model (kT's of 0.35, 1.01, and 6.0 
keV) does not fit the XMM-Newton PN data very well, 
especially in the energy regime below 0.7 keV, the low 
and high temperatures are almost identical. 

5.2. The X-ray light curve 

The lightcurve in Fig. |21 shows the development of the 
X-ray luminosity in the energy range 0.3 - 2.4 keV as de- 
termined from the 19 ROSAT observations, the Chandra 
observation from May 2000 and the new XMM-Newton 
data. The tendency over the first half year is character- 
ized by a decline with t^o.so (Zimmermann et al. 1993c, 
d), where t is the time since the outburst. Thereafter the 
lightcurve shows a bump and 5 years after the outburst 
the luminosity appears to return to the initial decline pro- 
file. The CHANDRA observation of 2000 and the XMM- 
Newton observation in 2001 are very close to the general 
^-0.30 decline. 

The mass loss rate of the wind of the progenitor can be 
expressed as M — iirpv^r'^, where v^^ is the wind veloc- 
ity. For a constant and homogeneous wind s has a value 
of 2. Under the assumptions that the density p of the 
circumstellar material at radius r is dominated by that 
wind and the shock velocity stays constant, the observed 
luminosity is roughly proportional to the square of the 
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Fig. 2. X-ray light curve showing the development of the 
X-ray luminosity in the 0.3-2.4 keV energy range. The 
dotted line indicates a decline rate with t~°-^°. The dashed 
lines result from model calculations involving time limited 
decrease and increase of the density in the ambient matter. 



density integrated over the emitting volume. In case of 
a constant and homogeneous wind the luminosity should 
show a time dependence oi oc t~^. The flatter time 
profile Qf ROSAT luminosity observed during 

the first half year, corresponds for a spherically symmet- 
ric scenario to a density profile with s=1.65. In a more 
detailed consideration in the framework of the formal- 
ism of the standard model the development of the X- 
ray luminosity produced in the forward shock is expressed 
by: oc T°-'^e^(3-2s){n-3)/{n-s) [(Fransson et al. 1996)| 

where r~" describes the radial dependence of the density 



of the SN ejecta. For n >> 10 it follows oc ^-(^''-s) 
and taking the observed slow decline with the value 

for s « 1.65 is confirmed, in good agreement with results 



published earlier by|(Immler et al. 2001) 



In the picture of the standard model the shock front, 
while running outward, is probing increasingly older peri- 
ods of the progenitor wind with a velocity roughly 1000 
times the wind velocity. Thus in 5 years ROSAT scanned 
about 5000 years of wind history. In our simple formalism 
the measured time dependence of the luminosity requires 
that the ratio 'mass loss rate/wind velocity' decreased dur- 
ing the last 5000 years in the life of the progenitor by more 
than a factor of 7. This scenario has been investigated in 
more detail by (Van Dyk 1994) using radio data and by 
(Immler et al. 2001) using the ROSAT data to calculate 
absolute mass loss rates. The latter authors also speculate 
that the change of the wind properties might indicate a 
transition from a red to a blue supergiant phase where 
mass loss rates are at least an order of magnitude lower 
and wind velocities may reach 1000 km/s. 

The bump in the X-ray light curve suggests a local in- 
crease in density above the general power-law profile of the 
circumstellar matter. It could be attributed to a change in 
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the wind parameters of the progenitor or to an asymmetry 
caused by a possible binary scenario. 

Interesting is a correlation with the expansion velocity 
of the SN shell that has been derived from high resolution 
radio images of SN1993J. In Fig. 6 of |(Bartel et al. 2002)| 
at about 350 days after the outburst the expansion rate 
of the radio image size shows a break and continues there- 
after at a slower rate. Between day 369 and 1655, the 
leading edge and the plateau of the bump, the expansion 
velocity decreases in the order of 28%. The total mass 
which gave rise to the observed X-ray luminosity bracket- 
ing the bump is estimated to 1 Mq. The total mass of the 
expanding shell covered up to day 350 was of the order of 
3 to 4 Mq. Both features, the radio deceleration break and 
the bump in the ROSAT light curve, independently indi- 
cate a flattening of the density profile in the circumstellar 
medium. Because the break in the expansion curve does 
not visibly affect the circularity of the radio images one 
may argue that the cause more likely can be attributed to 
a density change in the wind instead of an asymmetry by 
a binary scenario that would possibly affect the circularity 
of the images. 

The closer inspection of the X-ray light curve shows a 
rapid decline between day 225 and day 370, followed by 
a rapid increase until day 935. A similar behaviour could 
have occured between day 1126 and day 1464, but there 
is just one data point in that period, the upper limit of 
which is consistent with no change. The rapid decline af- 
ter day 1655 is very evident and recovery indicated by the 
Chandra data point at around day 2500 has occured. The 
existence of the latter two dips is statistically not too con- 
vincing, but they may be real as we show below. We sug- 
gest that each of the three dips indicates a density depres- 
sion which is followed by a rise of the density. If the shock 
wave runs into a sufficiently low density regime no more 
matter will be shocked and heated, the matter heated so 
far will simply expand. If the expansion is isothermal the 
light curve will go down with t^^; if the expansion is adi- 
abatic the decline will go like t^^/'^, as long as the tem- 
perature is higher than about 2 keV. This temperature 
is about what we would expect at around day 350 (cf. 
Fig. |3J). The falling dashed lines in figure 3 represent such 
expansion. After some time the shock wave encounters a 
jump in density to significantly higher values. The den- 
sity increase raises the emission measure and the X-ray 
luminosity. Initially the luminosity increases linearly with 
t when it hits a step like density jump. The rising lines 
in figure 3 represent such density jumps. The intersection 
between the rising and falling lines denote the position 
or epoch of the density jump. We have compared these 
epochs with those where the highest, discontinuous decel- 
eration of the shock wave has occured, which have been 
derived from the expansion rate of the radio images (cf. 
Fig. 6 of (Bartel et al. 2002) l. There are also just three 
such events in the radio images. The two sets of epochs 
agree surprisingly well with each other to within less than 
30 days. We conclude that these changes in the X-ray light 
curve reflect changes of the density in the ambient matter 
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Fig. 3. X-ray determined emission temperatures. The 
dashed lines outline the time evolution of the low and 
high temperature component using the standard model. 
For each component there are two lines bracketing the 
range of the temperatures allowed by the XMM-Newton 
data using the 2vmekal model. 



profile, which appears to show some repetitive pattern. 
Each of these density increases is preceded by some den- 
sity decrease like in a wave. Whether this characterizes 
the activity of the progenitor star as far as the mass loss 
is concerned remains to be seen. 

5.3. Spectral development 

We can now combine the XMM-Newton temperatures 
with the ROSAT, ASCA and GRO results. Let us assume 
for a moment that the initial high temperature, as mea- 
sured by GRO, originated from the forward shock, while 
half a year later the flux was dominated by emission from 
the reverse shock because of the higher electron densi- 
ties in that region. One can use the similarity solution of 
the standard model to express in general terms the ra- 
tio between the temperatures in the forward and the re- 
verse shocks that depends only on the parameters s and 
n,Tr./Tf = (3 - s)V(n - 3)^ |(Fransson et al. 1996) 



With the wind density profile of s=1.65 derived from 
the light curve, and assuming that the 2 temperatures de- 
termined from the thermal 2-component fit to the XMM- 
Newton spectrum represent the temperatures in the for- 
ward and reverse shock regions, we obtain n=8.9 for the 
density profile of the ejecta. Taking into account the uncer- 
tainties of the temperatures as quoted for the PN 2vmekal 
model, then n results between 6.4 and 10.9. 

Entering s and n into the formula that describes 
the temporal behaviour of the shocks, Tr cc Tf oc 
^-[2(3- s)/(n-s)] (Fpansson et al. 1996) we can extrapolate 



the XMM-Newton temperatures backward in time and 
compare them with the ROSAT, ASCA and GRO results 
(see Fig. |2l . The temperature development of the forward 
shock component is roughly consistent with the tempera- 
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tures deduced from the GRO data, and also the low tem- 
perature development, matching the ROSAT data at day 
210, is consistent with its interpretation as coming from 
the reverse shock region. 

But it has to be repeated that, while the parame- 
ter s appears reasonably well determined from the light 
curve (at least for certain periods), the value of n is de- 
rived under the assumption of thermal equilibrium be- 
tween electrons and ions in the emission region. Estimates 
(see (Fransson et al. 1996) i, however, show that even near 
to the time of the explosion equipartition is not expected 
in the wind shock. 

One can not exclude that we see just emission from 
the forward shock. The hot material heated in the early 
days has cooled as it expands, both radiatively and adia- 
batically with the continuum temperature dropping faster 
than the higher ionization levels indicate, i.e. we have a 
delay of the recombination because of the long recombi- 
nation time scales involved. Thermal equilibrium by just 
Coulomb collisions and ionization equilibrium is no more 
reached at times later than about a few days after out- 
burst. Stellar wind matter heated at later stages as the 
forward shock proceeds is not in thermal equilibrium if 
not processes other than Coulomb collisions are at work. 
Therefore the electrons would stay relatively cool and also 
ionization equilibrium is far from being reached leaving 
the plasma underionized, which would explain the obser- 
vation of the low temperature component. 

With the sedov model we have tested such a scenario 
with emission from just one shocked region. But the rela- 
tively poor fitting results for the line complexes at higher 
energies show also that this model is not quite adequate. 

6. Conclusions 

The XMM-Newton spectrum is well fitted by a thermal 
2-component model with ionization equilibrium. The pre- 
dictions of the standard model describe reasonably well 
the observed development of the X-ray spectral temper- 
atures, involving both forward and reverse shocks. The 
equally acceptable fit with a vsedov model shows that at 
this stage one cannot exclude that only one shock compo- 
nent is at work. 

The strong interdependence between the temperature 
and the values for the abundances of the dominating ele- 
ments, as far as the high temperature component is con- 
cerned, makes it very difficult to interpret the abundance 
results. But there is clearly no need to involve high val- 
ues for the metal abundances, which would point to ejecta 
heated by the reverse shock. 

The new XMM-Newton data point in the X-ray 
lightcurve gives arguments to the assumption of a gen- 
eral decline rate of cx and it suggests that the 
bump in the light curve is probably due to an intermediate 
density fluctuation. The detailed inspection of the X-ray 
light curve shows that a density increase did not occur 
just once, but in some repetitive fashion. There are three 
dips in the X-ray light curve each one followed by an in- 



crease in luminosity. Each of the dips occurred very close 
to the times when the forward shock velocity underwent a 
dramatic deceleration as indicated by the expansion rate 
of the radio images. 

The observations show that spatial changes of the 
powerlaw density profile have no significant impact on 
the morphology of the radio images, that remain very 
much circular. This indicates that the bumps in the X-ray 
lightcurve are more likely caused by density changes in a 
spherically symmetric wind rather than by asymmetries 
introduced by a postulated binary companion. For most 
of the time we need a wind that in the past left more ma- 
terial per volume unit in the circumstellar environment 
than near to the time of explosion. 
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